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SYMBOLS 


Fr  =  Froude  number  defined  as  (V2  lgR0),  where  Fr»  1.  This  nondimensional 
number  is  the  ratio  of  the  inertial  force  to  the  gravity  force. 

G  =  current  collector  dimensional  radial  gap  (Fig.  1). 

H  =  current  collector  dimensional  axial  gap  (Fig.  1). 

L  =  dimensional  length  along  current  collector  rotor  at  free  surface  to  tip  of  rotor 
(Fig.  1). 

r*  =  dimensional  radial  coordinate  from  axis  of  rotor  where  origin  of  coordinate 
system  is  at  O'  (Fig.  1). 

x  =  ( r*-Ra)/H  =  nondimensional  radial  coordinate  from  free  surface  when  coor¬ 
dinate  system  is  located  at  point  O  (Fig.  1). 

z  =  nondimensional  axial  coordinate  where  origin  of  coordinate  system  is  at  O' 

(Fig.  1). 

a  _i _ a_ 

dr*  H  dx 

Ra  =  dimensional  radial  coordinate  from  axis  of  rotor  to  point  where  free  surface 

touches  rotor  where  origin  of  coordinate  system  is  at  O'  (Fig.  1). 

d  J _ d_ 

dz*  H  a  z 

P  =  nondimensional  pressure  normalized  by  qV2e  -  P*/qV2e. 

u  =  nondimensional  radial  velocity  component  normalized  by  V0c’/i  =  v* I V ot A . 

v  =  nondimensional  azimuthal  velocity  component  normalized  by  VQ  =  ve*!V0. 

V0  =  QR„,  where  Q  is  the  azimuthal  radian  angular  speed  of  rotor  (Fig.  1). 

w  =  nondimensional  axial  velocity  component  normalized  by  V0ih  =  v*!V0zh . 

6  -  azimuthal  angular  component  in  radians  (0  =  unit  vector).  Since  radians  are 

nondimensional,  there  is  no  superscript  (*). 

Um  =  dimensional  characteristic  of  meridional  velocities. 

Re  —  _  Reynolds  number  defined  in  this  work  which  is  considered  to  be 

M  an  0(1)  quantity. 

R*  =  dimensional  radius  of  curvature  at  free  surface. 

*  =  /oo(z)  +  */o2(z)  +  £2/o4(z)  +  •••  +  /v-,/,o(M  +  tV'/i^)  + 

tF~xfn[B,Z)  +  Fr~2f2n(6,z )  =  f(B,z)  =  nondimensional  perturbation 
expansion  of  free  surface. 
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SYMBOLS  (Continued) 


yoo (x,z)  +  tv^A)  +  £2yo4  +  ...  +  Fr  'y10Oc,e,z)  +  thFr  1  v,,(jc.©,e) 

+  tFr  vn(x,8,z)  +  Fr  y20(x ,8,z)  +  ...  =  nondimensional  perturbation 
expansion  series  for  generalized  variable  y. 

dimensional  force  term  in  Navier-Stokes  equations. 


nondimensional  surface  tension  term. 


n  =  [1  +  (Zoo')2]  1/2  [x-fooi]  -  unit  normal  along  free  surface  (Fig.  3). 

r *  =  R0  +  H  [/w  +  ^r~l  s>n  SF^J-'j  J  =  dimensional  free  surface 

to  0(1)  and  0(Fr~l). 


u  ~  uoo  +  F r  1  s*n  ®  U(x,z)  =  nondimensional,  nonaxisymmetric  asymptotic 
expansion  for  the  radial  velocity  component  perturbed  by  gravity  to 
0(1)  and  0(Fr"'). 

v  =  Vqo  +  F~]  sin  8  V(x,z)  =  nondimensional,  nonaxisymmetric  asymptotic 
expansion  for  the  azimuthal  velocity  component  perturbed  by  gravity  to 
0(1)  and  0{Fr-’). 

w  =  Woo  +  F~x  sin  8  W(x,z)  =  nondimensional,  nonaxisymmetric  asymptotic 
expansion  for  the  axial  velocity  component  perturbed  by  gravity  to  0(1)  and 
0{Fr~\ 


=  -P*6U  +  H 


=  dimensional  stress  tensor  in  an  incom¬ 


pressible  fluid  =  Ojj  =  qV0  e  Ojj,  where  oy  is  the  nondimensional  stress 
tensor. 


q  =  density  of  fluid. 

M  =  viscosity  of  fluid. 

g  =  acceleration  of  gravity. 

e  =  nondimensional  small  parameter  (i.e.,  e«1)  =  H/R0  used  in  perturbation 

expansions. 
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ABSTRACT 


Some  designs  of  liquid  metal  current  collectors  in  homopolar 
motors  and  generators  are  essentially  rotating  liquid  metal  fluids  in 
cylindrical  channels  with  free  surfaces  and  will,  at  critical  rotational 
speeds,  become  unstable.  An  investigation  at  David  Taylor  Research 
Center  (DTRC)  is  being  performed  to  understand  the  role  of  gravity  in 
modifying  this  ejection  instability.  Some  gravitational  effects  can  be 
theoretically  treated  by  perturbation  techniques  on  the  axisymmetric  base 
flow  of  the  liquid  metal.  This  leads  to  a  modification  of  previously 
calculated  critical  current  collector  ejection  values  neglecting  gravity 
effects.  The  purpose  of  this  report  is  to  document  the  derivation  of  the 
mathematical  model  which  determines  the  perturbation  of  the  liquid 
metal  base  flow  due  to  gravitational  effects.  Since  gravity  is  a  small 
force  compared  with  the  centrifugal  effects,  the  base  flow  solutions  can 
be  expanded  in  inverse  powers  of  the  Froude  number  and  modified  liquid 
flow  profiles  can  be  determined  as  a  function  of  the  azimuthal  angle. 

This  model  will  be  used  in  later  work  to  theoretically  study  the  effects  of 
gravity  on  the  ejection  point  of  the  current  collector.  A  rederivation  of 
the  hydrodynamic  instability  threshold  of  a  liquid  metal  current  collector 
without  gravity  effects  is  presented  in  an  appendix  to  this  report. 

ADMINISTRATIVE  INFORMATION 

This  work  was  a  cooperative  effort  between  DTRC  and  the  Department  of  Mechanical 
and  Industrial  Engineering  at  the  University  of  Illinois,  Urbana,  Ill.  61801.  The  work  was 
supported  by  the  DTRC  Independent  Research  Program,  Director  of  Naval  Research 
(OCNRIO),  and  administered  by  the  Research  Director  (DTRC0113)  under  Program  Ele¬ 
ment  61152N,  Project  Number  ZR00001,  Task  Area  ZR0230201,  Work  Unit  1-2712-125, 
Project  Title:  “Orientation  Effects  in  Liquid  Metal  Collectors.” 

INTRODUCTION 

The  development  of  superconductive,  direct  current  homopolar  motors  and  generators 
used  for  advanced  naval  machinery  offers  opportunities  to  achieve  significant  improvements 
in  naval  ships.  A  key  feature  in  the  successful  operation  of  these  advanced  machines  is  the 
confinement  of  the  liquid  metal  used  as  the  sliding  electric  contact  or  current  collector. 

When  a  fluid  with  a  free  surface,  e.g.,  a  liquid  metal  current  collector,  is  rotated  in  a 
cylindrical  channel  (Fig.  1),  a  critical  rotational  speed  is  reached  where  the  liquid  is  for¬ 
cibly  ejected  from  the  channel.  This  ejection  has  been  analyzed  and  correlated  with  the 
onset  of  the  Kelvin-Helmholtz  flow  instability  at  the  free  surface  resulting  from  a  protec¬ 
tive  cover  gas  flowing  over  the  liquid  metal.  The  linear  stability  analysis  performed  in 
previous  work  by  Woo  et  al. 1  consisted  of  two  sequential  parts:  first,  the  determination  or 
estimation  of  the  steady-state  flow,  which  is  the  base  solution  and  results  in  determining 
values  of  the  liquid  and  gas  azimuthal  velocity  VL  and  VG,  respectively,  at  the  free  surface 
and,  second,  the  determination  of  the  stability  at  the  free  surface.  In  the  theory  for  deter¬ 
mining  the  ejection  point,  the  inputs  are  the  values  of  VL  and  Vc  as  a  function  of  Q,  the 
rotational  velocity,  and  the  output  is  the  critical  value  of  Q  (i.e. ,  QCRIT)  at  which  ejection 
occurs  because  an  azimuthal  perturbed  wave  begins  to  grow  with  time.  This  analysis 
neglects  gravitational  effects  and  treats  only  perfectly  axisymmetric  base  flows  and 
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perturbations  so  that  one  azimuthal  position  is  entirely  equivalent  to  every  other  position. 
The  Appendix  presents  a  clearer,  more  precise,  and  somewhat  different  derivation  of  the 
hydrodynamic  instability  threshold  without  gravitational  effects  of  a  liquid  metal  current 
collector  than  those  presented  by  Woo  et  al. 1  However,  the  previous  results  were  confirm¬ 
ed  by  the  new  derivation. 

The  objective  of  this  investigation  is  to  determine  the  role  that  gravity  has  on  modify¬ 
ing  the  ejection  stability  at  different  orientations  of  the  current  collector.  This  information 
will  allow  the  design  of  liquid  metal  collector  systems  which  will  operate  in  a  variety  of 
orientations  and,  thus,  different  gravity  loadings.  The  overall  gravity  investigation  being 
performed  at  DTRC  is  divided  into  experimental  and  theoretical  parts.  An  experimental  ap¬ 
paratus,  whose  geometry  is  similar  to  a  current  collector  and  whose  orientation  with 
respect  to  the  direction  of  gravity  can  be  changed,  was  built  and  experiments  were  per¬ 
formed  to  determine  the  effects  of  gravity.  This  work  will  be  correlated  with  the  present 
effort  in  a  later  report.  The  other  part  of  the  investigation  was  to  determine  the  effects  of 
gravity  on  the  previously  developed  theoretical  stability  analysis  by  Woo  et  al}  Gravita¬ 
tional  effects  enter  into  two  parts  of  the  stability  analysis.  In  the  first  part  of  the  previous 
theory,  gravity  modifies  the  base  flow  solutions  by  making  the  flow  nonaxisymmetric.  The 
liquid  is  accelerated  or  retarded  at  various  azimuthal  positions  0,  thus  affecting  the  critical 
ejection  angular  velocity  Q  and  defining  a  preferred  angle  0  for  ejection.  In  the  second 
part,  gravity  adds  another  force  to  the  perturbed  wave  which  destabilizes  at  the  top  of  the 
collector  but  which  stabilizes  at  the  bottom  of  the  collector.  The  objective  of  this  report  is 
to  document  the  derivation  of  the  mathematical  model  to  determine  the  effect  of  gravity  on 
the  axisymmetric  current  collector  flow  profile.  This  mathematical  model  will  be  used  to 
develop  the  theory  for  the  critical  ejection  rotational  velocity  with  gravity  in  later  works. 
The  model  is  hydrodynamic  and,  thus,  does  not  include  the  magnetic  field.  This  effect 
will  be  investigated  in  later  work.  Since  gravity  is  a  small  force  compared  with  centrifugal 
effects,  the  base  solutions  can  be  expanded  in  inverse  powers  of  the  Froude  number  and 
modified  liquid  profiles  can  be  determined  as  a  function  of  0.  The  surface  tension  effects 
are  included  in  the  boundary  conditions  to  0(1). 

Other  sources  on  current  collector  magnetohydrodynamic  flow  theory  are  available.2  6 

CURRENT  COLLECTOR  CONFIGURATION 
DISCUSSION  OF  PROBLEM 

Figure  1  shows  the  geometry  of  a  typical  open-channel  current  collector  in  a  gravita¬ 
tional  field.  The  current  collector  basically  consists  of  a  disk  rotating  in  a  cylindrically 
shaped  groove  filled  with  a  small  amount  of  liquid  metal  to  provide  an  electric  contact 
between  the  disk  and  the  groove.  The  liquid  metal  is  protected  from  chemical  attack  by  an 
inert  cover  gas,  which  is  also  in  motion.  The  rotor  has  a  radius  Ra  +  L  and  rotates  in  the 
azimuthal  direction  around  the  z  axis  with  a  constant  angular  velocity  Q.  The  angle  0  is 
measured  from  a  horizontal  radius  in  the  direction  of  rotation  (see  Fig.  2).  The  axial  gap 
width  is  designated  as  H,  the  radial  gap  as  G.  In  the  current  collector  system,  the  rela¬ 
tively  small  dimensions  H,  G,  and  L,  with  respect  to  R0  have  the  same  order  of  magnitude 

designated  as  H  =  O  ( H )  =  0(G)  =  O  (L).  The  aspect  ratio  ^  will  be  a 

nondimensional  parameter,  much  less  than  one,  that  will  be  of  considerable  interest  in  the 

0  V Jit 

problem.  A  Reynolds  number  shall  be  defined  as  Re  =  _ 2- — ,  which  is  considered  to 
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be  of  0(1).  Here,  the  equation  V0  =  QRa  is  the  velocity  of  the  rotor  at  the  free  surface, 
which  q  and  ^  are  the  density  and  viscosity  of  the  liquid  metal. 

For  steady-state  axisymmetric  flow,  all  derivatives  with  respect  to  azimuthal  angle  9 
and  time  t  are  zero;  however,  in  this  problem  the  gravitational  field  perturbs  the  steady- 
state  axisymmetric  flow,  and  derivatives  with  respect  to  6  are  retained. 


DIMENSIONAL  HYDRODYNAMIC  EQUATIONS  WITH 
GRAVITATIONAL  PERTURBATION 

The  dimensional  hydrodynamic  equations  for  fully  developed  nonaxisymmetric  flow  in 
a  gravitational  field  can  be  expressed  as 
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and 
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(3) 


(4) 


Equations  1  through  3  are  the  components  of  the  Navier-Stokes  equations  in  cylindrical 
coordinates  with  a  gravitational  perturbation.  Equation  4  is  the  equation  for  incompressible 
flow.  Figure  2  shows  the  gravitational  field  in  the  cylindrical  coordinate  system  (r*,0).  9 
is  measured  from  a  horizontal  radius  in  the  direction  of  rotation. 
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The  dimensional  variables  designated  by  a  superscript  (*)  are  defined  as  follows 

v*  =  dimensional  radial  velocity  component 

v  g  =  dimensional  azimuthal  velocity  component 

v*  =  dimensional  axial  velocity  component 

r*  =  dimensional  radial  coordinate  (r  =  unit  vector) 

8  =  azimuthal  angular  component  in  radians  (0  =  unit  vector).  Since  radians  are 

nondimensional,  there  is  no  superscript. 
z*  =  dimensional  axial  coordinate  (z  =  unit  vector) 

P*  -  dimensional  pressure 

q  =  density  of  fluid 

p  =  viscosity  of  fluid 

g  =  acceleration  of  gravity. 
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Fig.  2.  Gravitational  force  in  r*  and  9  components. 


In  examining  Fig.  1,  it  is  noted  that  the  liquid  metal  is  contained  on  the  outer  rim  of 
the  rotating  disk.  Therefore,  a  local  radial  coordinate  Hx  =  r*-Ra  can  be  defined  at  the 
point  O,  where  the  interface  between  the  liquid  and  gas  touches  the  rotating  disk.  From 
this  definition  of  x  the  following  partial  differential  relationship  holds: 


3 _  =  Id 

dr*  Hdx 


(5) 


CHARACTERISTIC  QUANTITIES  AND  NONDIMENSIONALIZATION 

The  azimuthal  velocity  ve*  must  vary  from  zero  at  the  stator  to  VG  =  Qr  at  the  rotor. 

Taking  Ra  to  be  a  characteristic  radius,  the  characteristic  azimuthal  velocity  is  VQ  =  QRa. 

The  meridional  flow,  (/.  e. ,  secondary  flow),  with  the  radial  velocity  v*  and  the  axial 

'  *2 

velocity  v?*,  is  driven  by  the  centrifugal  force  due  to  the  azimuthal  motion  gVg  r  .  The 

r  ~  Dv* 

centrifugal  force  accelerates  the  meridional  circulation  until  either  the  inertial  term  q _ 1_ 

I"  v  *  "I  d* 

or  the  viscous  term  p  V*2vr*  -  _jL_  balances  the  centrifugal  force.  In  this  problem  there 

L  r*“  J 

is  interest  in  Reynolds  numbers  much  greater  than  one  (/?e»l).  Since  the  Reynolds 
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number  is  defined  as  the  ratio  of  the  inertial  force  to  the  viscous  force,  the  inertial  term  is 

(Dv  *  ov  *2\ 

q _ 1_  =  ~JL_  )  .  This  balance  determines 

Dt*  r*  ) 


the  characteristic  velocity  Um*  for  the  meridional  velocities  vr*  and  v*.  The  following 
characteristic  quantities  can  be  determined: 

d? 


L  \h) 


Dt* 


dr* 


and 


v0*  ~  V  =  QR  —  ~  _L 

VB  o  XIXo  »  * 

r*  Rn 


Thus,  the  force  balance  is  the  following  characteristic  quantity: 


QU 


*  2 


H 


where 


U  *  = 


qk; 

K 

JL 


V  =  z  V 

r  o  c  r  O  • 


(6) 


(7) 


(8) 


(9) 


a  p 


Balancing  the  pressure  gradient  -  - against  the  inertial  term  one  obtains  qV0  z  as  the 

dr 

characteristic  pressure  difference  due  to  the  centrifugal  force. 

Equations  1  through  4  can  be  made  nondimensional  by  using  the  following  nondimen- 
sional  variables. 


vr*  =  vorhu 


= 


v*  =  vyw , 

P*  =  qV2zP  , 


2*  =  Hz  ,  ~  =  — 


dz* 


H  dz 


r*  =  R0  +  Hx  =  z~xH{\+zx) 

d  ^  i  a 

dr*  H  dx 

The  following  equations  have  the  following  nondimensional  form  for  small  values  of  z. 


du  .  i a  du  3 u 

u —  +  (I+ex)  lz*v - 1-  w - (I+ejc)  v 

dx  30  dz 


-1..2 


(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 


=  -  ?JL  -  (gll  \  Sin  e  +  Re 

V  v} 


dx 


+  t  ( 1  4-tr ) 


1 T— 

Idx2 


_i  du 

+  E(1+£X)  - 

dx 


a  2u 
ae2 


a2u 

dz2 


e2(  1  +zx)~2u—2z  3/2(  1  +zx)~2 


3v1 

aej  ’ 
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(18) 


3  v 


U—  +£,/5(l+£jr)  1  V— -  +  wj—  +  e(l+ex) 
ax  do  dz 


3  v 


3v 


-l 


uv 


=  -  ii,2(l+ix)~l  —  -  l^-f 
dB  ' 


cos  0+Re~l  \—2 

\Vo)  ldx 


__  i  d  V  j  _ n 

+  f(l+£X)  —  +  £2(1+£.x)  2  2 

3jc  302 


dzv  3  v 
dz2 


-  e2(l  +ex)~2v 


+  2£5/2(1+£jc)-2 


du 

30  J  ’ 


3  w  w  ,  3  w  3  w  3  P 

“77  +  £  (1+“>  v  77  +  “77  =  '  77 


(19) 


and 


+  Re 


du 


'[ 


32w 
3*  2 


+  £(l+£Jt) 


3  w 


dx 


+  £2(1+£j:) 


_•>  32w 


30" 


+ 


32w 

dz2 


]■ 


—  +£(1+£X)  'u  +  ^(l  +  EJt) 

dx 

where 

u  =  u(x,B,z)  , 


-l 


3  v 
~30~ 


3vv 

3z 


=  0  , 


(20) 


v  =  v(x,0,z)  , 
w  =  w(x,6,z )  . 

The  Reynolds  number  was  defined  as  Re  = 


QUm*H  _qV0Hi* 


to  be  of  0(1)  in  this  problem.  The  Froude  number  is  defined  as  Fr  = 


and  is  considered 
,  the  inertial 


sK 


force  divided  by  the  gravity  force.  In  this  problem,  the  Froude  number  is  much  greater 
than  one  (Fr » 1). 


BOUNDARY  CONDITIONS 

At  all  solid  surfaces  of  the  current  collector  stator  and  rotor,  the  nondimensional 
secondary  flow  components  are  zero. 

u  =  u(x,6,z)  =  0  (at  solid  surfaces),  (21) 

and 

w  =  w(jc,0,z)  =  0  (at  solid  surfaces).  (22) 

At  the  stator  surface,  the  nondimensional  primary  flow  component  is  zero. 

Thus, 

v  =  v(x,6,z )  =  0,  at  stator,  (23) 

and  at  the  rotor  surface  the  nondimensional  primary  flow  component  is 

v  =  v(x,9,z)  =  ( 1  -hex),  at  rotor.  (24) 
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In  the  problem  described  in  this  report,  the  nondimensional  velocity  components  are  func¬ 
tions  of  x,8,  and  z  not  x  and  z,  as  in  the  axisymmetric  flow  case  considered  earlier  by 
Audet  et  al.6  This  is  because  the  gravitational  field  perturbation  makes  the  flow  in  the 
work  discussed  in  this  report  nonaxisymmetric. 

The  boundary  conditions  at  the  free  surface  must  be  specified.  The  surface  oc  in  the 
current  collector  configuration  in  Fig.  1  is  a  surface  of  separation  between  two  immiscible 
fluids:  the  inert  cover  gas  and  the  liquid  metal.  The  first  condition  to  be  satisfied  at  the 
free  surface  requires  that  the  velocities  of  the  fluids  in  the  direction  of  the  unit  normal  at 
the  free  surface  be  equal  to  each  other  and  equal  to  zero.  The  second  condition  requires 
the  forces  which  fluids  exert  on  each  other  at  the  free  surface  be  equal  and  opposite. 

These  conditions  are  written  respectively  as 

vini  =  Y  •  n  =  0  (25) 

and 

Ojjtij  (liquid)  +  a^rij  (inert  gas)  =  0.  (26) 

The  unit  vector  for  the  fluid  is  along  the  outward  normal  vector.  The  unit  normal  vec¬ 
tor  for  the  gas  is  in  the  opposite  direction.  In  order  to  simplify  the  boundary  conditions 
surface  tension  effects,  which  are  important,  have  been  neglected  at  this  point  of  the  work 
but  will  be  considered  later  in  this  report.  Neglecting  stresses  in  the  inert  cover  gas  at  the 
free  surface  equation  (26)  simplifies  to 


OijTlj  =  0, 

(27) 

and  * 

*  “ "  6,1  +  “  ( »v  +  ■«£-)  ■ 

(28) 

The  stress  tensor  components  o will  be  made  nondimensional  later  in  this  report. 


PERTURBATION  EXPANSIONS  OF  FUNDAMENTAL  VARIABLES 

The  coupled  boundary  value  partial  differential  equations  and  boundary  values 
discussed  previously  have  two  small  fundamental  parameters,  t  and  F~ 1 ,  for  perturbation 
expansions.  An  estimate  of  these  parameters  for  a  typical  current  collector  configuration  is 
shown  in  Fig.  1 .  If  the  axial  distance  is  H  =  3  x  10-4  meters  and  the  distance  along  the 
rotor  surface  from  the  z  axis  to  the  liquid  interface  is  Ra  =  0.3  meter,  the  parameter  i  = 
HIR0  =  10~3  and  ~  0.0316  =  3.16X  10~2.  For  a  typical  speed  of  3, (XX)  revolutions 
per  minute  or  314.2  radians  per  second,  the  velocity  Q 0R0  =  94.25  meters  per  second. 

The  Froude  number  Fr  =  V02/gR0  is  calculated  to  be  3,021.3,  Fr-1  equals  3.31  x  10~4 
and  F~l  «1.  The  Froude  number  is  a  nondimensional  parameter.  The  Reynolds  number 

Re  =  is  considered  as  an  0(1)  parameter  and  does  not  enter  the  asymptotic 

M  Vi  -i 

analysis.  In  an  asymptotic  expansion,  z  and  Fr  are  considered  to  approach  zero. 
ASYMPTOTIC  EXPANSION 

The  limits  as  zh  and  F~ 1  -*■ 0  are  considered  independently,  so  a  double  asymptotic 
expansion  for  the  velocity  v  is  required.  Thus, 

Y  =  loo  +  £*Yoi  +  £Y02  +  £3/2Y03  +  •  •  •  (29) 
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+  fT'vio  +  +  *Fr-lYn  +  £3V.v13  +  ■  •  ■ 


+  Fr  2Y20  +  e*Fr  ^ ¥21  +  &r  2  *22  +  •  •  • 
+  Fr  3  Y30  +  ^Fr~3Y 31  +  •  •  • 


+  Fr  4i’40  +  •  •  • 

+  .  .  .  . 

This  explicit  form  applies  for  each  component  of  velocity  u,  v,  and  w  and  for  the  pressure  P. 

Gravitational  effects  first  enter  the  expansion  in  Eq.  29  on  the  second  line.  Each  term 
in  the  first  line  v0j,  for/  =  0,  1,  2,  . . .  is  independent  of  gravitational  effects. 

Without  gravitational  effects,  the  flow  is  perfectly  axisymmetric.  Therefore,  v00,  v0I , 
Vq2 ,  y03,  etc.,  are  all  independent  of  0.  In  addition,  the  axisymmetric  version  of  the 
governing  equations  does  not  involve  th,  only  t.  For  the  first  row  of  this  double  series, 
the  only  small  parameter  is  e,  not  tA,  so  that  v01  =  0,  v03  =  0,  v05  =  0,  etc. 

Therefore, 

V  =  Vqo  (X,z)  +  £v0 2  (x,z)  +  c2vM  (x,z)  +  . . .  (30) 

+  Fr~lYio  (x,e,z)  +  e“Fr-1v11  (x,0,z)  +  tF~ 1  y ,2  (x,0,z) 

+  Fr~2Y20  (x,d,z)  +  ...  . 


Only  two  terms  in  this  series,  Vqq  (jc,z)  and  y,0  (x,6,z),  where  Vqq  (x,z)  is  the  ax¬ 
isymmetric  base  solution  will  be  computed.  The  base  solution  neglects  gravitational  effects 
and  assumes  £«1,  while  v10  (x,0,z)  is  the  first  perturbation  of  the  base  solution  due  to 
gravitational  effects,  v^,  will  be  ignored  even  though  e  is  comparable  to  F~l.  Mathema¬ 
tically,  this  double  series  is  valid,  as  e-»0  and  independently.  Therefore,  the 

expressions  for  v10  and  can  be  independently  derived.  Each  term  in  the  series  depends 
only  on  terms  above  and  to  the  left  of  it.  Physically,  ev^  (x,z)  may  be  a  larger  perturba¬ 
tion  than  Fr~lvl0  (x,6,z),  but  the  former  is  still  axisymmetric,  while  the  latter  represents 
the  first  deviation  from  axisymmetiy.  The  location  of  the  ejection  point  is  determined  by 
the  deviation  from  axisymmetry,  regardless  of  how  small,  because  this  deviation  deter¬ 
mines  which  point  becomes  unstable  slightly  before  another  point. 

Substituting  u  =  ( x,z ),  v  =  Vqq  (x,z),  w  =  (jc,z),  and  P  =  (x,z)  Eqs.  17 

through  19  and  the  boundary  conditions  result  in  the  following  system  of  coupled  base 
equations: 


du 


‘oo 


‘oo 


oo 


3jc 
d  v 


+  w 


du 


00 


00 


00 


dx 


+  W 


dz 

dv 


00 


00 


dz 


2  _  dF00 
-  v°o  -  ~  a7~ 

+  Re~l 

[a2«oo  | 

La*2 

*V] 

3*2  J  ’ 

(31) 

-  Re~x  P2"00  + 

a2vool 

(32) 

Lax2 

3Z2  J 
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dw, 


* 00 


oo 


dx 


,  dwoo  _  3Poo 
+  woo  - - -  - - +  Re 


-l 


and 


du 


oo 


+ 


3  w. 


00 


dz 


=  0 


dz 


[ 


d2wpo 

dx2 


^  2>v00~| 
dz2  J  ’ 


(33) 


dx  dz 

This  is  the  same  result  obtained  by  Audet  et  al.  ,6  except  they  made  the  pressure  p*  non- 


(34) 


dimensional  by  P  = 


H 


-,  which  differs  from  the  factor  which  makes  the  pressure 


nondimensional,  Qv2e,  (i.e.  P*  =  qv2lP),  by  a  factor  of  the  Reynolds  number.  The  work 
of  Audet  et  al.6  is  fully  axisymmetric  and  does  not  include  the  effects  of  gravity. 

BOUNDARY  CONDITIONS 

The  governing  nonlinear  partial  differential  equations  determine  (x,z),  (x,z), 

w(K)  (jc,z),  and  P^  (x,z)  for  specified  boundary  conditions.  The  following  boundary  condi¬ 
tions  are  specified  at  the  current  collector  surfaces.  At  all  solid  surfaces  both  rotor  and 
stator  are 

Woo  (x,z)  =  Woo  (x,z)  =  0  .  (35) 

v’oo  (x,z)  =  0  at  the  stator  surfaces,  and  Vqq  (x.z)  =  1  at  all  rotor  surfaces.  For  the  boun¬ 
dary  conditions  at  the  free  surfaces,  the  components  of  the  stress  tensor  are  required  and 
are  normalized  by  qV2e..  With  Eqs.  17  through  20  made  nondimensional  and  the  coor¬ 
dinate  system  being  changed  from  r,0,  and  z  to  x,8,  and  z 


*  _ 


=  qK,  «»,•, 


(36) 


The  components  of  the  dimensional  stress  tensor  oi}*  from  dimensional  cylindrical  coor¬ 


dinates  to  nondimensional  coordinates  o,,  are: 


=  -  P*  +  2p 


3v/ 

dr* 


_i  3“ 

or  =  -  P  +  2Re  1  — 
dx 


}ee 


=  -/>*  +  2p  — 


1  dvg* 


de* 


* 


(37) 

(38) 


°ee  *  -  p  +  2Re~lt/‘(l+rx) 


o*  =  -P*  +  2M 


dv^ 

dz* 


3w 

=  -  P  +  2Re  1  —  , 
dz 


dv 

W 


+  tu 


(39) 


°r%  =  °6r*  =  M  (  — 


1  d  v. 


30 


+ 


w 

dr* 


Jx8 


=  °6x  =  Re 


-1 


[■-  (S 


+  e^(l  +£jc)' 


(40) 


t-0]- 


10 


a  =  Q 
zr  rz 


/£v  +  »v\ 

\3r*  3z*  ^ 


(41) 


°*z  ~  °« 


and 


_,  T  3k  3w  "I 

=  /te  1  -  +  -  , 

|_  3z  dx  J 


*  _ 


'zfl  =  °9z  =  M 


W 

dz* 


3v/ 

33 


(42) 


?ez 


=  °ze  =  ^ 


i  \c*  dv 
dz 


u  _|  3  M> 

+  ‘  (1+“>  -ST 


] 


In  order  to  derive  the  stresses  on  the  free  surface,  the  free  surface  location  must  be 
known.  However,  the  free  surface  x  =  /  (0,2)  is  another  unknown  of  the  problem. 
Therefore,  the  free  surface  will  be  expanded  in  a  perturbation  series.  Thus, 


*  =  /oo(z)  +  */<»(*)  +  £ 7o4(z)  + 

+  Fr~lf\o(8>z)  +  eAF-lfn(9,z)  +  tf~lFn(6,z)  +  ... 
+  Fr~2f2o(0>z)  +  •  •  •  =  Mz)  ■ 


(43) 


P  (x,0,z),  u  (. x,6,z ),  v  ( x,9,z ),  and  w  (. x,B,z )  are  expanded  in  the  same  form  as  x  =  J(9,z). 
All  of  these  unknowns  are  governed  by  the  partial  differential  equations  (Eqs.  17  through 
20),  the  boundary  conditions  at  the  rotor  and  stator  surfaces  and  the  four  boundary  condi¬ 
tions  at  the  free  surface  (normal  velocity  and  three  components  of  surface  traction  equal 
zero).  The  solution  for  each  unknown  depends  on  the  parameters  t,  F~l,  Re,  and 
geometry.  Each  variable  can  be  replaced  by  its  double  asymptotic  expansion  in  t  and  F~\ 
where  the  coefficient  functions  still  depend  on  Re  and  the  geometry  but  are  independent  of 
e  and  F~x.  Like  the  other  variables,  / (0,2)  has  a  first  row  which  is  independent  of  0 
(axisymmetric)  and  involves  only  t,  not  e  . 

The  boundary  conditions  on  u^,  Vqq,  Wqq,  and  at  x  =  f^iz)  are: 

nji  +  nzw  =  0,  (44) 

and 


for  i  =  (x,6,z), 
where  the  nondimensional  surface  tension  term 


(45) 


(t)~SC-(£) 


and  R* 


equals  the  dimensional  radius  of  curvature  along  the  free  surface.  Rc*  is  considered 
positive  when  the  surface  is  convex  in  the  gas  phase.  Here  y  is  the  surface  tension  and  a 
=  VlQVftH  is  a  dimensionless  surface  tension  coefficient.  The  surface  tension  term  is 
necessary  for  numerical  stability  for  numerical  calculations  to  zero  order.  The  dimensional 
radius  of  curvature  is  expressed  as 

11  -/of' 


HR„ 


U  +  ifm’n 


21  3/2 


The  unit  normal  along  the  free  surface  n  is  defined  as 

n  =  (l+(/oo')2r*  li-fo'o  z]  , 


(46) 


(47) 
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where 
and 


nx  =  [1+  O^)2]-*,  «9  =  0 


and 


*z  =  ~/oo  n+  obo)2]^. 
x  =  radial  unit  vector, 
z  —  axial  unit  vector, 


f,  _  dfoo(z) 
/o°  ~Jz 


Therefore,  the  base  boundary  solutions  at  the  free  surface  are  at  x  =  fw(z) 
are  as  follows: 

“oo  =/oo  woo  » 

9voo  _  f ,  3voo 
dx  00  dz 


~  Re  pm  -  Re 


and 


foo  RePoo  +  foo  Re 


(  «  \  ,  „  a“oo  _  ,,  / 3 “oo  ,  9woo\ 

[X  j  sr) 

(«r) 


+  9“oo  +  dwoo  -  2 f '  dw°° 
dz  dx  100  dz 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


MODIFICATION  AND  DISCUSSION  OF  BOUNDARY  CONDITIONS 

Traditionally,  two  linear  combinations  of  the  last  two  equations  are  taken.  The  first 
combination  eliminates  the  pressure  *oo: 


fo'o(z)  [Eq.  59]  +  [Eq.  55]. 

(57) 

Th“s>  o  a  o  o 

9 “00  L  3u00  L  3w00  _  ,s,^2  3u00 

2/00  dx  +  dz  +  dx  ffoo)  3* 

(58) 

+  (Too)2  l”™  *  *"'»  . 

dx  dz 

At  x  =  fwiz),  the  normal  component  of  the  velocity  vn  equals  zero.  v„  can  be  differen¬ 
tiated  along  the  surface  using  the  chain  rule,  and  the  result  can  be  used  to  simplify  the 
condition  (Eq.  58)  further.  TTie  condition  on  v„  is 

Uoo(foo(z),z)  «  fo'o(z)  Wootfooiz),  z).  (59) 


Differentiating  with  respect  to  z,  the  following  equation  is  obtained: 


foo 


d“ooC foQ’Z) 
dx 


+  3uQo(foQyZ) 

dz 


~  foo  wOo(foO’Z) 


(60) 


+  (foo)2 


3woo(foo>z) 

dx 


+  foo 


3wQo(foQ’Z) 

dz 
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Subtracting  Eq.  60  from  Eq.  58,  the  following  boundary  condition  is  obtained: 


+  foo  woo  =  0  -  at  *  =  /«,(*)  . 

Equation  61  represents  the  stress  condition,  o„,  =  0,  where  (n,6,t)  are  orthogonal  cur¬ 
vilinear  coordinates;  the  unit  normal  f  is  tangent  to  the  surface  x  =  f^iz)  in  each  0  = 
constant  plane,  as  shown  in  Fig.  3.  If  the  surface  is  a  straight  line  in  each  0  =  constant 
plane,  then 


z 


Fig.  3.  Free  surface  x  =  /qqIz)  with  normal  vector  n  and  tangential  vector  t. 
dv 

Since  vn  =  0  here, -  =  0,  which  is  equivalent  to  subtracting  Eq.  60  from  Eq.  58.  The 

dt 

boundary  condition  for  a  straight  line  surface  is 


dv.  fa  ,  a  1  ,,  (63) 

— - U  +  (foo)  ]  I  — —  foo  ~ —  (foouoo  +  woo)> 

dn  \_dx  dz  J 

which  reduces  to  the  first  four  terms  in  Eq.  61  for  a  straight  line  surface  with  =  cons¬ 
tant  and  Jqq"  =  0.  The  fifth  term  in  Eq.  61,  f<x)"woo>  represents  the  effect  of  surface  cur¬ 
vature  in  each  0  =  constant  plane.  At  each  point  on  the  free  surface,  polar  coordinates  are 
taken  in  a  0  =  constant  plane  with  the  origin  at  the  local  center  of  curvature.  In  these 
polar  coordinates,  (R,a),  the  surface  is  at  R  =  Rc  =  local  radius  of  curvature,  and  a 
changes  along  the  surface.  Now  the  equation  becomes 


°R„  =  Re 


1  dvR  , 

Rc  da  3R 


(64) 
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where 


1  dvR  _  dvn  __ 

R  da  dt 
and  the  final  new  term,  - 


3v„  _  3v, 

dR  dn  ’ 
v 

—  ,  represents  the  extra  term  due  to  curvature.  Since  um 

Rc 


/oo'  Wqq  at  the  surface,  the  following  equation  is  obtained: 


(65) 


V„  =  V,  =  u  +  (/o(,)2r‘A  (foo  «oo  +  vvoo)  (66) 

=  [1  +  (/oo)2]"^  . 

The  local  radius  of  curvature  is  provided  by 

=  -  foo  I1  +  </oo)2r3/2,  (67) 

so  that  the  extra  term  due  to  surface  curvature  is 


“  ~  =  /oo  [1+  (fo'ofV^oo-  (68) 

" c 

The  sum  of  the  straight  line  surface  condition  (Eq.  63)  and  the  extra  term  due  to  surface 
curvature  (Eq.  68)  provides  the  boundary  condition  of  equation  (61). 

The  other  linear  combinations  of  the  two  equations  involving  Pm  are  obtained  from 
Eqs.  55  and  56: 

fob  (Eq-  56)  -  (Eq-  55). 

Thus, 

p00  =  (  -  2Re~l  n  +  (/oo)2]'1  f—22  -  fob  d-~  (69) 

\  Rc )  t  dx  3  z 

-foo  ’at  x=fco(z)- 

This  equation  represents  the  boundary  condition  normal  to  the  surface 


and 


vn  —  U  +  (/oo)2]  W  [Moo  foo  woo]  •  (72) 

Curvature  terms  do  not  enter  into  this  equation.  In  the  following  work,  Eqs.  61  and  70 
shall  be  used  for  the  boundary  conditions. 

The  boundary  value  problem  for  (x,z),  uw  (x,z),  (x,z),  (x,z),  and/*,  (z) 

is  the  nonlinear  problem  currently  being  solved  numerically  by  Audet  et  al. ,  for  the  DTRC 
current  collector  program  for  various  high  values  of  the  Reynolds  number.  Audet  et  al.6 
have  developed  a  numerical  scheme  for  the  solution  of  the  zeroth  order  Navier-Stokes 
equations  in  order  to  calculate  the  equilibrium  flow  distribution  in  the  liquid  metal  region 


14 


of  a  current  collector  (Fig.  1).  The  computer  program  calculates  problems  for  Reynolds 
numbers  of  up  to  1 ,000.  The  solution  of  the  finite  difference  equations  is  based  on  a 
Newton-Raphson  method.  A  mathematical  method  is  being  developed  to  calculate  the  free 
surface  shape  from  the  boundary  conditions.  Audet  et  al.6  are  in  the  process  of  calculating  a 
variety  of  velocity  profiles  for  a  number  of  current  collector  geometries;  henceforth,  it  can 
be  assumed  that  this  base  solution  is  known  and  the  first  perturbation  can  be  considered 
due  to  gravitational  effects,  i.e.,  P10  ( x,8,z ),  «j0  ( x,6,z ),  v,0  ( x,9,z ),  w10  (x,0,z),  and/10  (0,z), 


GRAVITATIONAL  PERTURBATION  EXPANSION 


M  a  ,  a2 

In  Eqs.  17  through  20,  the  0  derivatives  appear  as  z  - or  z  — =- 

30  dGr 


The  first  term 


with  0  variation  is  Fr-1vi0,  so  that  the  0  derivative  terms  are  all  z'hF~x  or  z2F~l,  which 
are  negligible  compared  to  the  F~l  term.  Therefore,  0  enters  through  the  gravitational 
terms,  but  no  0  derivatives  enter  so  that  the  solution  can  be  obtained  in  each  0  =  constant 
plane  independently  of  all  other  0  =  constant  planes. 

The  two  gravitational  terms  are  -F~l  sin  0  in  the  radial  equation  and  -F~x  z'h  cos  0 
in  the  azimuthal  equation.  The  former  drives  the  perturbation;  the  latter  term  is  negligible 
in  this  problem. 


Therefore,  the  first  gravitational  perturbation  terms  can  be  written  as: 


m10  (x,8,z)  =  sin  6  U  (x,z)  ,  (73) 

v10  (x,9,z)  =  sin  0  V  ( x,z )  ,  (74) 

w10  (r,0,z)  =  sin  0  W  ( x,z )  ,  (75) 

/io  (9>z)  =  sin  0  F  (z)  .  (76) 

The  effects  of  gravity  are  maximum  and  opposite  at  the  top  (0  =  90°)  and  bottom  (0  = 
270°);  they  are  negligible  at  0  =  0°  and  0  =  180°.  When  the  asymptotic  expansions 

v  =  Vqo  (x,z)  +  F~x  sin  0  V  {x,z)  (77) 


are  introduced  into  Eqs.  17  through  20  and  the  terms  0(e),  0  ( F~lz'k ),  and  0  ( F~ 2)  are 
neglected  and  the  equations  governing  Uqq,  Vqq,  Wqq,  and  P0 0  are  subtracted  and  the  re¬ 
mainder  is  multiplied  by  Fr  (esc  0),  then  a  set  of  linear  equations  is  obtained  governing 
U,V,W,  and  P,  which  have  variable  coefficients  involving  the  known  solutions  u^,  v^, 
and  woo. 


Thus, 


3  U  d  //qo 

"“a7  +  JTU+W‘ 


,  duoo 
°°dz  dz 


(78) 


»  dP  ,  n  f  d2U  d2U~ 

dx  [dx2  dz2] 


The  - 1  in  Eq.  78  is  derived  from  -  sin  0  and  is  the  only  inhomogeneous, 

driving  term  in  the  boundary  value  problem;  the  known  functions  come  first: 
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(79) 


where 


(80) 


(81) 


U  =  U  (x,z)  , 

V  =  V{x,z)  , 
and 

W  =  W  (. x,z )  . 

The  following  boundary  conditions  apply  to  Eqs.  79  through  81.  Equation  82  is  true  at  all 
solid  surfaces,  both  stator  and  rotor. 


Thus, 


U  =  0,  V  =  0,  and  W  =  0.  (82) 

At  the  free  surface,  the  following  overall  conditions  remain: 


v,n i  =  0,  where  i  =  x,  Q,  z  , 


(83) 


and 


°ijni  = 


(84) 


The  four  boundary  conditions  expressed  in  Eqs.  83  and  84  will  be  evaluated  in  the  next 
section  to  0  (F,~r)  to  show  the  effects  of  gravity.  The  effects  of  surface  tension  are 
neglected  at  0  (F~ 1 )  because  the  surface  tension  effects  are  not  needed  for  numerical 
stability  as  in  the  case  of  the  zeroth  order. 


GRAVITATIONAL  PERTURBATION  BOUNDARY  CONDITIONS 
UNIT  NORMAL  TO  FREE  SURFACE 


The  dimensional  coordinate  at  the  free  surface  r*  is  defined  to  0(1)  and  0  (Fr  '). 


Thus,  . 

r*  =  Ra  +  H  T/qo  +  F~x  sin  6  F  (-j- 

)]' 

00 

.  L  V  /  V 

where 

x  =  /oo(z)  +  F~x  sin  9  F(z)  , 

(86) 

and 

r*  =  R0  +  Hx  ,  zH  =  z*. 

(87) 
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The  function  |  can  be  defined  as 


t  <r\e,z-l  =  r*-R0-H  L,  (0  +  F~'  Sto  8 F  (01 

_  r  _  r\  ai_ »  £ _ _  _ _ i  .. ‘*u  a  /r  "K 


(88) 


where  E  =  0  at  the  free  surface  perturbed  with  0  (Fr  )  gravity  terms.  From  the  definition 
of  the  gradient  vector  V,  the  normal  to  the  free  surface  is 


n  |  =  ££_ 

4=0  dr * 


-  1  34 

r  +  —  — 
r*  ae 

4=0 


9  + 


4=0 

1  17  “1 


3£_ 

3z* 


z 

4=0 


=  r  -  £  (l+tx)-1  F,_1  cos  0  F(z)  6 

-  Ifooiz)  +  Fr~l  sin  0  F'(z)]  Z  • 
The  unit  normal  is  then  expressed  as 


n  =  V£ 


/|?«|  I 

=0/  1  4=0 


=  n 


,o)  (l?{Uo 


-1 


where  the  components  nx>  ng,  and  nz  are: 


nx=  Yt 


'4=0 


-1 


ne  =  -t(l+ex)-lF-llVt 


-i 


cos  0  F(z)  , 


4=0 


and 


nz=  -  V*| 


I  \  f  > 

|  Too 

4=0 


/oo(z)  +  Fr  1  sin  9  F 


r,(z)J  . 


(89) 


(90) 


(91) 

(92) 


(93) 


DERIVED  GRAVITY  BOUNDARY  CONDITIONS 

The  boundary  condition  Eq.  83  will  be  evaluated  and  can  be  expressed  as 

unx  +  vng  +  wnz  =  0  (94) 

at  the  free  surface.  Substituting  nx,  ng,  and  nz  from  Eqs.  91  through  93  and  keeping  terms 
to  0(1)  and  0(Fr-1),  results  in  the  expression 

u  =  [fooOO  +  Fr~l  sin  0  F'(z)]w  , 

at  X  =  /oo(z)  +  F~l  sin  9  F(z)  ,  (96) 

substituting  u  =  ( x,z )  +  F~l  sin  9  U  ( x,z )  and  w  =  Wqq  ( x,z )  +  Fr~l  sin  9  W  (x,z) 

and  using  Eq.  96  for  x  at  the  free  surface  provides 

“oo  Ooo(z)  +  Fr~ 1  sin  0  F(z),z)  (97) 

+  Fr~i  sin  0  U  Ooo  (z)  +  Fr-1  sin  0  F  (z),  z) 
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=  I/o6<*)  +  Fr~X  sin  0/6o  (2)] 

X  {w00  (/oo(Z)  +  Fr~l  sin  0  F(z),  Z> 

+  F,.-1  sin  0  W  (/oofe)  +  Fr_1  sin  9  F(z),  z)}  . 

Careful  examination  of  Eq.  97  reveals  that  it  is  not  an  asymptotic  expansion.  A  true 
asymptotic  expansion  must  consist  of  powers  or  functions  of  the  small  parameter  or 
parameters  ( F~‘)  times  coefficient  functions  which  are  independent  of  F~x .  Terms  such 
as  Uoof/oofz)  +  F~l  sin  9  F(z),  z)  depend  on  F~l. 

Therefore,  the  Taylor  series  expansion  for  u ^  is  introduced,  which  is 


«00  (/oo(z)  +  Fr  1  sin  0  Ffe),  Z)  =  «00  (/oo(z),  Z) 

+  F,-'  sin  «  m  —  -,X--  Z- 

dx 

+  'AF~2  sin20  [F(z))2  + 


(98) 


dx1 

Also,  analogous  series  for  U,  and  W  are  developed.  These  Taylor  series  expansions 


are  substituted  in  Eq.  97,  and  terms  for  0  ( Fr  )  are  neglected  to  obtain: 
“oo  (/oo(z)»  z)  +  F~l  sin  9  |V(z)  z) 


(99) 


[Ft 


),  z ) 


+  V  (/oo(z) 

'(zHo(/oo(z),  z)  +  fooiz)F(z) 
+  /oo(z)^(/oo(z) 


/oo(z)  Woo  (/oo(z),  z)  +  Fr  1  sin  0 
aw00(/00(z),  z) 


0* 


h  z)  . 


Subtracting  Eq.  53,  the  0(1)  part  of  the  boundary  condition,  from  Eq.  97  and  multiplying 
by  Fr  esc 0  results  in  the  gravity  boundary  condition  at  0  ( F~l).  Thus, 


a “oo  (/oo(z).  z) 
dx 


F (z)  +  0(/oo(z),  z) 


(100) 


=  Woo  (/oo(z),  z)  F'(z)  +  /ri  (z)  a->-(^°(z)’  Z)  F(Z) 


0a: 


+  /oo(z)H'(/oo(z),  z)  ■ 


The  three  boundary  conditions  from  Eq.  84  (o,yi^=0)  will  be  evaluated  at  0(  Fr  '). 
These  three  conditions  are: 

+  «e°98  +  nz°ez  *  0  * 

+  «9°x8  1  " 


and 


»*°M  +  «8°z8  +  nz°zz  =  0  • 


■«-0, 

(101) 

.  =  0 , 

(102) 

=  0  . 

(103) 
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These  three  expressions  are  analogously  compared  to  Eq.  94  at  0(  Fr  ')  to  show  the  ef¬ 
fects  of  gravity.  They  are,  respectively: 


x  a  Voof/oofz),  z)  ,  dV(foo(.z),  z) 
r*(z)  — : -  -t- 


0JC2 


(104) 


dx 


=  foo(z)  F'(z) 


dvoo(foo(z),  z) 
dx 


+  /oo(z)  F (z) 


d2voo(foo(z)>  z) 
dxdz 


+  foo(z)  v(foo(z)>  z) 

+  ffe)  £> , 

3z 


_  iteF(z)  3^oo(-0°a)’-?)  -  RePifooiz),  z) 
dx 

+  2  F(z)  d2Uoo(foo(z),  z)  +  2  dUjfwiz),  z)  = 
djc2  3jc 


+  /oift)  JTO 


aMoo(/oo(z),  z) 
3* 


+  /oo(2)  TO 


a2Moo(/oo(^)» z) 
3.r3z 


dz 


+  /oo(^)  TO 


a  ^oo ( Too (z ) ;  z) 
Sjc2 


+  foo(z)  U(foo(z),  Z) 

+  F'(z)  au°°(/oo(z)’  z) 


+  f>(-)  3h,oo(/oo^’  z) 
dz 


P,(z)  a«oo(/oo (Z),  z)  +  F(z)  a2Moo(/oo (Z),  z) 
dx  dxdz 

+  dUjfooiz),  z) 
dz 


(105) 


(106) 


19 


DISCUSSION  AND  CONCLUSIONS 

The  base  axisymmetric  flow  solutions  of  the  current  collector  axisymmetric  flow  field 
are  designated  as  Uqq,  vm,  and  w^,,  which  are  the  nondimensional  radial,  azimuthal,  and 
axial  velocity  components,  respectively.  Gravity  pertuibs  the  axisymmetric  base  flow  solu¬ 
tions  of  the  current  collector  flow  field  by  making  the  flow  nonaxisymmetric.  Gravity 
accelerates  or  retards  the  base  flow  at  various  azimuthal  positions.  The  nondimensional, 
nonaxisymmetric  asymptotic  expansions  for  the  velocity  components  perturbed  by  gravity 
are  determined  to  be  of  the  form: 

M  =  Moo  +  ?r~'  S‘n®  >  (107) 

V  =  Vqo  +  F~'  sine  V(x,z)  , 
and 

w  =  Woo  +  l  s*n®  W(x,z)  , 

where  the  gravity  perturbation  term  is  of  order  0(Fr~’)  and  where  Fr  is  the  nondimensional 
Froude  number  (the  ratio  of  the  inertial  force  divided  by  the  gravity  force).  The  gravity 
term  is  also  a  function  of  sin0,  where  0  =  0°  is  along  the  horizontal.  The  functions  U,V, 
and  W  are  only  functions  of  x  and  z,  and  not  0. 

Work  is  presently  being  performed  by  Audet  et  al.6  to  determine  the  base  flow  cur¬ 
rent  collector  solutions  without  gravity  perturbations  Mqo,  Vqq,  and  Wqq  numerically  from 
Eqs.  31  through  34  with  appropriate  boundary  conditions  including  the  surface  tension 
effects  (Eqs.  53  through  56).  The  first  condition  at  the  free  surface  requires  that  the 
velocities  of  the  fluid  vectors  in  the  direction  of  the  unit  normal  at  the  free  surface  be 
equal  to  each  other  and  to  zero.  The  second  condition  requires  the  forces  which  the  fluids 
exert  on  each  other  to  be  opposite.  In  this  work  the  stresses  in  the  inert  cover  gas  are 
assumed  to  be  small  enough  to  be  neglected.  The  functions  U,  V,  and  W,  the  functions  for 
the  gravity  perturbations,  are  governed  by  four  coupled  linear  partial  differential  equations 
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(Eqs.  78  through  81),  which  are  functions  of  the  base  solutions  Uqq,  Vqq,  and  %  as 
variable  coefficients.  The  gravity-perturbed  boundary  conditions  are  determined  by  Eq.  82 
and  the  linear  boundary  conditions  (Eqs.  100  and  104  through  106).  This  complex  system 
of  equations  will  be  used  in  future  work  to  calculate  the  effect  of  gravity  on  the  base  flow 
of  the  current  collector  in  order  to  determine  the  effects  of  gravity  on  the  linear  stability 
analysis  that  determines  the  collector  ejection  speed.  Woo  et  al.1  have  performed  the  linear 
stability  analysis  without  gravity. 
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APPENDIX 

REDERIVATION  OF  LINEAR  STABILITY  ANALYSIS  WITHOUT 

GRAVITY  EFFECTS 


INTRODUCTION 

This  Appendix  presents  a  clearer,  more  precise,  and  somewhat  different  derivation  of 
the  hydrodynamic  instability  threshold  of  a  liquid  metal  current  collector  than  that 
presented  by  Woo  et  al. 1  This  derivation  also  confirms  the  theory  developed  previously  for 
the  hydrodynamic  ejection  threshold.  At  high  rotor  speeds,  the  liquid  metal  current  collec¬ 
tor  is  subject  to  flow  instabilities  leading  to  almost  instantaneous  ejection  of  the  liquid 
metal  from  the  current-carrying  region  and  disrupting  the  operation  of  the  device.  The 
basic  geometry  considered  in  this  study1  is  the  classic  tongue-and-groove  illustrated 
schematically  in  the  cross  sectional  view  in  Fig.  4,  in  which  an  inner  disk  is  fitted  into  the 
circumferential  groove  of  the  outer  disk.  Generally,  the  grooved  piece  is  the  stator,  and  the 
disk  is  the  rotor.  In  Fig.  4,  the  radial  gap  is  of  width  d;  the  axial  gap  of  length  a  between 
the  disk  and  the  grooved  piece  are  exaggerated  for  clarity.  The  volume  in  the  radial  gap  is 
filled  with  liquid  metal,  which  is  set  into  rotational  motion  by  the  drag  from  the  rotation 
of  the  inner  disk  and  which  is  maintained  in  the  grooved  volume  by  centrifugal  force. 

Thus,  an  electrical  current  path  is  maintained  between  the  disk  (rotor)  and  the  groove 
(stator).  To  minimize  both  electrical  and  viscous  power  losses,  the  volume  of  the  liquid 
metal  occupies  a  minimum  of  the  axial  gap.  As  a  result,  when  the  liquid  metal  is  fully 
disturbed  in  the  collector,  there  is  a  free  surface  in  the  axial  gap  between  the  liquid  metal 
and  a  nonreactive  cover  gas.  Generally,  the  radius  to  the  liquid  metal  surface  R  is  much 
greater  than  d  and  a.  The  cover  gas1  is  also  in  rotational  motion  due  to  drag,  and,  as  with 
the  liquid  metal,  the  azimuthal  flow  is  a  function  of  the  position  along  the  axis  because  the 
boundary  conditions  require  the  cover  gas  flow  to  match  both  the  moving  and  stationary 
surfaces.  The  centrifugal  forces  cause  secondary  circulation  flows  in  both  the  liquid  metal 
and  the  cover  gas.  The  meridional  secondary  flows  form  vortices  superposed  on  the 
azimuthal  primary  flow.  In  general,  the  stability  of  the  current  collector  flow  field  in 
response  to  perturbations  is  of  concern.  The  stability  of  the  liquid  metal  collector  from 
hydrodynamic  first  principles1  will  be  discussed  in  this  Appendix.  Magnetohydrodynamic 
effects  on  the  current  collector  ejection  will  be  considered  in  later  work  by  the  authors. 
Further  information  can  be  obtained  from  the  efforts  of  Eriksson,7  who  pioneered  earlier 
work  on  this  problem.  The  variables  used  in  the  Appendix  are  defined  as  needed  for  the 
reader,  and  are  not  necessarily  those  used  in  the  text. 


PRIMARY  FLOW  FIELDS 


The  dynamics  of  the  current  collector  system  in  Fig.  4  are  described  by  the  dimen¬ 
sional  Navier-Stokes  equation  and  the  continuity  equation  in  cyclindrical  coordinates. 


(108) 


(110) 
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STATOR  LIQUID  METAL 


Fig.  4.  Typical  current  collector  configuration.1 


(HD 


(112) 

For  steady  state,  axisymmetric  flow,  derivatives  with  respect  to  time  t  and  azimuthal 
angle  are  zero.  The  fluid  has  a  viscosity  \i  and  a  density  q. 

The  dimensional  variables  are  defined  as: 

(«,vfl,w)  =  components  of  the  velocity  field  in  the  radial,  azimuthal,  and  axial  directions.; 
(r,0,z)  =  components  of  cylindrical  coordinates  in  radial,  azimuthal,  and  axial  directions. 

P  =  scalar  pressure  component. 

The  conducting  fluid  flowing  between  the  rotating  disk  and  the  stationary  groove  is 
contained  absolutely  except  at  the  interface  with  the  inert  cover  gas  between  the  azimuthally 
rotating  disk  and  stationary  sidewall.  Although  the  situation  will  become  unstable  for 
sufficiently  large  differencies  between  the  azimuthal  velocities  in  the  liquid  metal  and 
cover  gas,  stable  equilibrium  will  be  established  at  lower  azimuthal  velocities.  In  the 
mathematical  model,  the  radial  gap  between  the  stator  and  rotor  is  assumed  to  be  zero. 

The  gas  is  assumed  to  be  incompressible  and  occupies  the  region  0  <  r  <  R.  The  liquid 
metal  occupies  the  region  R  <  r  <  R  +  d.  The  gas  and  liquid  base  velocities  are  assumed 
to  be  similar  to  two  different  rigid  body  rotations  with  no  axial  variations.  VL  and  VG 


where  the  substantial  derivative  DIDt  and  the  Laplacian  V2  are 


£  -  —  +  u  — 
Dt  dt  dr 


+  w- 


r  30 

and  the  continuity  equation  is 


dz 


•,V2  = 


1 

+  - 


3_  +  _l_ 


dr2  r  dr  r2  dQ2  dz‘ 


1  d(ru )  1  3  w 

- 1 — _  + - i  +  =  0  . 

r  dr  r  30  dz 
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represent  the  interface  liquid,  and  gas  primary  flow  velocities  at  r  =  R,  respectively, 
which  are  constants.  The  primary  solutions  for  the  velocity  and  pressure  for  the  liquid 
i  =  L  and  gas  i  =  G,  respectively,  are: 


V6L  ~ 


V9G  -  Vc 


+  T/R 


(113) 

(114) 


where  T  is  the  surface  tension  between  liquid  and  gas. 

The  primary  pressure  in  the  liquid  PL  and  gas  PG  is  obtained  from  the  radial  com¬ 
ponent  of  the  Navier-Stokes  equation  (Eq.  108)  by  equating  the  centrifugal  force  term 
2 

to  the  radial  gradient  of  the  pressure,  neglecting  the  other  terms,  and  integrating 
between  0  and  r.  After  substituting  V6i  =  V{  ^ ,  integration  provides 


i  vy 

—  P(r)  =  -I-  constant. 

q  2 R2 


(115) 


Arbitrarily  setting  the  pressure  of  the  liquid  PL  to  0  at  interface  r  =  R  determines  the  con¬ 
stants  in  both  equations  for  the  liquid  and  gas  and  results  in  Eqs.  113  and  114. 


PERTURBATION  OF  VELOCITY  AND  PRESSURE 

The  perturbation  expansions  are  restricted  to  have  no  axial  variation,  — —  =  o, 


and  no  axial  velocity,  vz  =  0.  The  expansions  for  the  liquid  velocity  vL,  liquid  pressure 
PL,  gas  pressure  vG,  and  gas  pressure  PG  are  the  two  first  orders  in  e. 

Thus, 


V*,  =  t  vrL(r,e,t)  r  + 

I  V1  (— ^  +  £  v9 L  ( r,e,t ) 
l  \R  J  J 

§  - 

(116) 

Q'Vr2  /  r2  ^ 

2  (*-') 

|  +  «  PL  ( r,B,t )  , 

(117) 

vG  =  e  'VcM.O  l  + 

[VG  +  £  ^c(r-0-O 

?  . 

(118) 

\  T 

)  +  —  +  e  Pc(r,Q,t )  , 

(119) 

where  t  is  a  small  parameter  much  less  than  1,  i.e.,  e«l. 

Since  gravity  is  neglected,  the  problem  has  axisymmetric  symmetry.  The  centrifugal 
force  in  the  system  due  to  the  rotational  motion  of  the  fluid  is  the  dominant  equilibrium 
force.  Viscous  terms  are  neglected  in  this  analysis  because  this  is  an  inviscid  analysis. 

Only  0(e)  terms  will  be  kept  in  this  analysis,  and  nonlinear  0(e2)  terms  and  those  of  higher 
order  will  be  dropped  because  this  is  a  linear  analysis  for  small  perturbations  for  £«1. 
The  hydrodynamic  stability  of  the  collector  interface  at  r  =  R  between  the  liquid  and 
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cover  gas  in  relative  motion  was  studied  by  the  perturbed  Navier-Stokes  equation  and  the 
perturbed  continuity  equation  to  first  order  in  t  for  both  fluids. 

Thus, 


a  vr,  v.  a  v  -  2V:  1 

a>, 

—  +  -1  —  — L  h, 

_  i 

t 

(120) 

a t  R  30  R  J 

dr 

1 

aPt 

(121) 

=  —  — 

> 

dt  R  R  dd  J 

r 

dd 

vri  i  a  Vo. 

+  —  + - -  =  0  . 

(122) 

r  r  dB 

and 

a  Ki 

dr 

There  is  no  axial  component  to  the  Navier-Stokes  equation  since  there  is  no  axial 
velocity  perturbation.  These  equations  were  derived  by  substituting  Eqs.  116  through  119 
in  Eqs.  108  through  112  and  retaining  the  terms  to  order  t. 

BOUNDARY  CONDITIONS  AT  NONLINEAR  FREE  SURFACE 

The  free  surface  position  for  the  primary  flow  at  r  =  R  is  perturbed  to 

r  =  R  +  tf(B,t)  .  (123) 

Before  the  boundary  conditions  for  the  linearized  case  are  considered,  the  fully 
nonlinearlized  boundary  conditions  will  be  derived.  The  nonlinear  free  surface  condition 
can  be  expressed  as 

r  =  F(9,t)  ,  (124) 

where  the  liquid  is  in  the  region 

r  >  F(0,t )  .  (125) 

The  velocity  of  the  liquid  vL  and  the  pressure  of  the  liquid  PL  are,  respectively: 


Y.i  =  vrL(r<e’0  l  +  vet(r,0,/)  9  * 


and 


pl  ~  . 

The  velocity  of  the  gas  vG  and  the  pressure  of  the  gas  Pc  are,  respectively: 


vc  =  vrG<>'’0’/)  i  +  vec(r>®>/)  9  > 


and 


PG  =  Pc(r’®’0  • 


(126) 

(127) 

(128) 
(129) 


The  stress  tension  o,;  in  an  incompressible  inviscid  fluid  takes  the  simple  form  for  the 
liquid  L  and  gas  G,  respectively, 


°'JL  =  ~  Pl  6ij  . 


and 


°UG  ~  PG  • 


(130) 

(131) 
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Since  there  are  no  shear  stresses,  there  are  no  restrictions  on  tangential  velocities  at  the 
free  surface.  Since  the  fluids  are  considered  to  be  inviscid,  the  liquid  metal  and  cover  gas 
can  slide  over  each  other  at  the  free  surface.  The  inviscid  nonlinear  boundary  conditions  at 
r  =  F(0,r)  are: 


ZL  •  1 


=  Yg 


=  V  n 


(132) 


and 

PG 

where 

n 


T 

PL  +  — -  (normal  stress  condition)  , 
Rc 

local  unit  normal  to  free  surface, 
velocity  of  free  surface  position,  i.e., 


dF(6,t)  r 


(133) 


(134) 


and 

Rc  =  local  radius  of  curvature  of  free  surface  defined  as  positive  when  radius  of 
curvature  is  on  the  gas  side  of  the  surface. 

The  local  unit  normal  to  the  free  surface  can  now  be  derived,  and  the  function  can  be 
defined  as 


<D(r,0,r)  =  r  -  F(6,t)  ,  (135) 

so  that  O  =  0  at  the  free  surface.  V0  is  perpendicular  to  the  free  surface  at  r  =  F(8,t). 
Therefore,  the  gradient  of  0  is 


1  dF(d,t) 

V0  =  r - 0  , 


r  30 


where 


V  =  ?—  +  0  —  — 


'  3r  ~  r  30 
The  gradient  of  F  at  r  =  F(8,t)  along  the  free  surface  is 


V0 


=  r 


1  dF(8,t)  g 

F(9,t)  30 


1  r  =  F(0,<) 

Thus,  the  local  unit  normal  along  the  free  surface  n  is  equal  to 


V0 

r  =  F(0,t) 

1  V0 1 

”1 

II 

r  - 


1 


dF(6,t )  q 


F(0,t)  30 


(136) 


(137) 


(138) 


(139) 


(140) 
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Using  Eq.  132,  the  nonlinear  velocity  boundary  condition  is  found  to  be 

vri  =  —  —  vdi  +  —  at  r  =  F(8,t)  .  (141) 

f  dd  dt 

The  radius  of  curvature  at  each  point  along  the  interface  r  =  F  (0, t)  must  now  be 
derived.  A  line  r  =  F  ( Q,t )  in  the  z  =  0  plane  is  now  considered  (Fig.  5).  Let  s  equal  the 
distance  along  the  curve  from  some  point,  s  —  0.  n(s)  is  the  local  unit  normal  at  each 
point  on  the  curve.  The  distance  ds  equals  Rcda,  where  a  is  the  angle  subtended  by  ds  . 


I  / 

// 

^  -« - CENTER  OF  CURVATURE 


Fig.  5.  Radius  of  curvature  at  each  point. 


28 


The  unit  normal  at  n(s  +  ds)  can  be  expressed  in  terms  of  n(s)  , 

Thus, 

.  -  dn(s)  i 

n(s+ds)  =  n(s)  +  ds  — - H  0  (ds)2  +  . . .  (142) 

ds 


in  a  Taylor  series  expansion.  By  vector  subtraction 
_  „  dn(s)  -) 

n(s+ds)  -  n(s)  =  — —  ds  +  0  (ds)2  +  ...  .  (143) 

ds 

Since  |  n(s)  |  =  1,  n(s)  can  change  direction,  i.e.,  rotation,  but  not  length  (Fig.  5). 

Therefore,  is  perpendicular  to  n(s).  From  similar  triangles  the  identity 
ds 


dn(s) 

ds 


ds  =  |  n(s)  |  da 


(144) 


can  be  obtained,  from  which  the  radius  of  curvature  Rc(s)  can  be  derived 


1 

Rc(s) 

since 


dn(s) 

ds 


(145) 


dn(s) 

ds 


(146) 


The  curve  in  the  z  =  0  plane  is  r  =  F  (8,t).  The  two-dimensional  metric  in  cylin¬ 
drical  coordinates  (r,0)  is 


ds  =  [(dr)2  +  r2(de)2]'A  at  r  =  F(9,t)  . 


(147) 


dF 

Since  dr  =  —  dd  at  a  particular  time,  ds  can  be  expressed  again  as 
oa 


The  derivative  — 
ds 


is,  simply, 


and  the  derivative 


dn  . 
—  is 
ds 


-'h 


(148) 


(149) 


dn  _  dn  dd  _  2  +  /  9F\  2  'h  dn 

ds  ~  dB  ds  L  \30 /  J  de 


(150) 


Actually,  the  n  derived  from  Eqs.  139  and  140  is  a  function  of  8  and  r,  but  the  curvature 
at  each  point  on  the  curve  for  one  particular  time  is  being  addressed,  so  t  is  fixed  here. 
Thus, 
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[f2  +  (lf)1 "  "  ("tf)  -] 

derivative  of  n  with  respect  to  9  is 

[>  +  /— VT  f—  ?+rtl 

de  L  \  de )  J  La0'  'J 


(151) 


(152) 


x 


where  it  is  remembered  so  that 

dr  »  ,  d6 

—  =  9  ,  and  —  =  -  r  . 
d9  '  d9 


(153) 


From  Eqs.  150  and  152,  the  radius  of  curvature  of  the  curve  at  r  =  F  (9,t)  at  a 
given  instant  is 


1 


dn 

ds 


\  30  )  J  1  de 


(154) 


■[-©ri'-i'-©! 


Now,  the  only  question  is  whether  the  correct  sign  for  this  convention  is  obtained.  If 

F(6,t)  -  R,  Eq.  154  provides  which  is  correct.  Therefore,  the  boundary 

Rc  R 

condition  for  the  liquid  and  gas  pressure  at  the  interface  is 

PG  =  PL+  T/Rc  at  r  =  F(6,t )  .  (155) 

These  are  the  boundary  conditions  for  a  nonlinear  analysis. 

LINEAR  BOUNDARY  CONDITIONS  AT  A  FREE  SURFACE 

The  nonlinear  free  surface  was  represented  as  r  =  F(0,r).  A  linear  representation  of 
the  free  surface  would  be  r  =  R  +  The  nonlinear  free  surface  velocity  boundary 

condition  was  derived  previously  as 
1 


3  F  ,  3  F 

—  vfll-  +  — 

ae  e‘  dt 


at  r  =  F(6,t), 


where  i  =  L  or  G 


and  v,  =  tvri  ( r,9,t )  r  +  Vi  ^  +  ev6i  (r,0,r)  8. 


The  boundary  condition  to  order  0(e)can  be  expressed  as 

i  df(9,t ) 

£vri  [R  +  zf(9,t),e,t]  =  [R  +  zf(6,t)] 

OU 

x  jV,  |\  +  +  tv9L  [R  +  e/(0,f)]j 


(156) 

(157) 

(158) 


+  i 


1  + 

dt 


30 


(159) 


Introducing  the  well  known  Taylor  series: 
vri  [R  +  tf(8,t),9,t ]  =  vri  (R,8,t) 

+  t  ^  (R,e,t)  +  Vi  £2/2  («.M  +  •  •  • . 

3r  dr2 

and 


vei  [R  +  ef(9,t),9,t]  =  vgi  ( R,9,t ) 

+  ef(8,t)  ^  (. R,8,t )  +  Vi  e2/2(0,O  ^4^  (*,9,0  +  •  ■  •  . 
dr  dr 2 

Substituting  these  expressions  in  Eq.  158  and  neglecting  terms  of  O  (e  2): 


K  df(8,t ) 


vri  ( R,e,t )  =  -i 

r  dd 

for  i  =  L  or  G. 


dt 


(160) 


(161) 


It  is  of  note  that  the  boundary  condition  is  applied  to  vri  at  r  =  R  because  the  Taylor 
series  was  about  the  point  r  =  R. 


The  expression  for  the  radius  of  curvature  to  order  t  is  obtained  from  Eq.  154  and  is 
equal  to 


(e2)  , 


i  i  e  r d2/(9,t )  ,  „  l 

Rc  R  R?  I  +  (  °J  +  0 

where  the  following  expressions  were  used: 

[f2  +  (“i^) 2]  =  ^  +  2JW>')  +  °<e2)  - 

[f2+(lF)T”  =  T-Fm')  +  0(£j)' 


3  F\  2 


(f) 


-  F 


d2F 

dd2 


„  3  7(0,0  „  2 

=  -  tR  -  -  +  0(e2)  . 

dd2 


(162) 


(163) 

(164) 

(165) 


The  linearized  normal  stress  condition  at  the  free  surface  is  obtained  from  the 
nonlinear  surface  condition 


where 

Pl  = 

and 

PG  = 


+ 


T 

R 


(166) 

(167) 

(168) 


After  simple  substitution  and  the  use  of  the  linearized  free  surface  condition  r  =  R  +  tf(9,t), 
the  expression  to  order  e  for  the  reciprocal  of  the  curvature,  (Eq.  162)  and  the  Taylor 
series  expansion  of  the  linearized  normal  stress  condition  at  the  free  surface  is 
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(169) 


Cg^g 


QlYl 


mo  +  pc  (Ms)  = 

R  R 


mo 


t  r  3  /(0.o  i 

+  PL(R,8,t)  -  —  ,  +  mo  . 

p2  Lae2  J 


PERIODIC  SOLUTION  IN  THE  AZIMUTHAL  ANGLE  AND  TIME 

Since  0  is  the  azimuthal  angular  coordinate,  the  following  substitution  shall  be  made: 
P,(r,8,t)  =  Pj(r)  exp  ( irnid  +  zcur)  .  (170) 

All  integer  values  shall  be  considered  and  the  value  of  n  corresponding  to  the  first  unstable 
mode  shall  be  found.  The  n  for  the  first  unstable  mode  is  very  large,  thus,  n  can  be 
treated  as  a  continuous  variable  rather  than  a  discrete  variable  with  only  integer  values. 
Therefore,  the  discrete  variable  nn  shall  be  replaced  with  the  continuous  variable  kR. 
Explained  differently,  the  first  unstable  mode  has  a  wavelength  2n Ik  for  the  first  unstable 
mode  k.  The  circumference  of  the  interface  is  2nR.  If  the  wavelength  is  much  smaller  than 
the  circumference 

—  «2nR  -  kR  »1  ,  (171) 

k 

then  the  disturbance  is  unaffected  by  the  finite  circumference  and  is  the  same  as  a  distur¬ 
bance  on  an  infinite  length  surface. 

The  simple  transformation  is  introduced: 

r  =  R  +  x  and  ~ —  =  —  ,  (172) 

dr  dx 

where  no  assumptions  on  the  size  of  x  have  been  determined  yet  in  the  work.  The  per¬ 
turbed  quantities  are  assumed  to  be  of  the  exponential  form 


vn(r,0,r)  =  Uri(x)  exp  [i(kR6  +  cur)]  , 

(173) 

v0l(r,0,r)  =  Uei(x)  exp  [i(kR9  +  cur)]  , 

(174) 

P;(r,0,r)  =  P,(x)  exp  [i(kR8  -1-  cur)]  , 

(175) 

/(0,r)  =  A  exp  [i(kR8  -1-  cur)]  , 

(176) 

where  i  =  L  or  G  (liquid  or  gas) . 

Substituting  Eqs.  173  through  176  the  perturbed  Navier-Stokes  equations  to  order  t 
results  in  the  following  set  of  ordinary  coupled  differential  equations: 


and 


?,  [«•  +  kvt)  Uri(x)  -  l/0t(*)j  =  - 

?,  J^t(w  +  kVt)  UBi(x)  +  -  Uri(x)  j  = 


dx 

ikR 
C R+x ) 


dUJx)  +  J^x)_  +  JkR_  ^  =  Q 


dx 


(R+x)  (R+x) 


(177) 

(178) 

(179) 
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IV  2V 

In  Eq.  177  the  term  ikViUri  »  — -  U6i,  and  in  Eq.  178  the  term  ikVtUei  »  — ■  Uri. 

R  R 

In  each  case,  the  smaller  term  is  neglected.  These  inequalities  were  determined  by  noting 
kR»\. 

The  radial  penetration  of  the  disturbance  from  the  free  surface  shall  be  shown  to  be 

2n  ^ 

comparable  to  the  azimuthal  wavelength,  i.e.,  — ,  —  is  comparable  to  k.  Noting 

k  dx 

that  R  »x,  the  term  Uri/R  is  neglected  in  Eq.  179.  Equations  177  through  179  are  now 
simplified  to  the  following  set  of  coupled  differential  equations: 


and 


/Ci(co  +  k  Vt)  Uri  = 

_ dPi 

dx 

(180) 

iQi(&  +  k  Vt)  U6i  = 

-  ikPj  , 

(181) 

dU 

— -  +  ik  U6i  =  0  . 

(182) 

dx 


From  these  three  equations,  the  following  simple  relationships  are  obtained: 
k  dx 

D  -  'C‘  (co  +  k  Vt)  dUri 


and 


d2U, 


dx‘ 


r 

ft  =  k2U, 


dx 


(183) 

(184) 

(185) 


The  general  solution  to  the  ordinary  differential  equation  (Eq.  185)  is 

Uri  =  Bxexp  ( kx )  +  B2  exp  (- kx )  ,  (186) 

where  it  is  assumed  that  k  >  0.  For  the  gas  (i  =  G),  the  following  relationships  hold: 

-  R  <  x  <  0  and  — —  =  k  -*■  bx  =  ,  where  Ax  equals  the  radial  extent  of  the  distur- 

dx  k 

bance  from  the  free  surface.  In  the  gas  RI(Lx)  =  kR  »1  and  x  is  located  effectively  in 
the  region  -  00  <jc  <0  and  UrG~*0  as  x  =  -  °°,  implying  Bz  =  0.  Therefore,  using 
Eqs.  173  through  176  and  Eqs.  183  and  184,  the  following  relationships  can  be  obtained 
for  the  gas  phase: 


UrG  (x)  =  Bx  exp  (kx)  , 
uec  (x)  =  iB\  exP  (**)  . 

PG  (x)  =  -  Bx  (co  +  kVG)  exp  (kx)  . 
k 


(187) 

(188) 

(189) 


For  the  liquid  (i  =  L),  the  following  relationships  hold  0<x< 

2n 

distance  d  »  Ax  azimuthal  wavelength  —  is  assumed,  then  kd  »  1 . 


d,  —  ~  kd.  If  the  radial 
Ax 

Therefore,  d  »  Ax  and  x 
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is  effectively  in  the  region  0  <x<  °°  and  UrL  0  as  x  -*■  impling  fl,  =  0.  Therefore, 
from  Eqs.  183  through  185,  the  following  equations  are  obtained  from  the  liquid  phase: 

UrL  (-*)  =  B2  exp  ( -kx )  ,  (190) 

UeL  (*)  =  ~  exP  (_£*)  >  (191) 

and 

Pr  (jc)  =  —  (co  +  kVL)  B2  exp  (~kx)  .  (192) 

k 

The  solutions  have  been  obtained,  except  for  a  constant  in  Uri  Ox),  Uei  (x),  and  P,  (x), 
where  /  =  L  or  G  for  the  gas  and  liquid  phases  in  the  current  collector  system. 

Using  the  linearized  velocity  condition  at  the  free  surface  for  the  liquid  L  and  the  gas 
G  and  after  appropriate  substitution,  the  following  boundary  conditions  are  obtained: 

UrG  (0)  =  fl,  =  i  (to  +  kVG)  A  ,  (193) 

UrL  (0)  =  B2  =  i  (o>  +  kVL)  A  ,  (194) 


where 

-  =  —  k2  R2  Aexp  [j  (kRQ  +  co  r)]  ,  (195) 

3d2 


and  /  (0,f)  in  the  surface  term  is  neglected. 


The  linearized  pressure  condition  at  the  free  surface  (Eq  169)  has  the  following  form 
after  appropriate  substitution: 

Qlvl  2  Gg^g 


~(w  +  k  VG)2  +  ^(co  +  kVL )2  = 
k  k 


+  k 2T , 


(196) 


where 

co 


(Qg^g  +  QlYl) 

(Ql  +  GC) 

*  W  ,  (v„  - 

(Gg  +  Gt) 


(197) 


Equation  197  was  obtained  from  Eq.  196  by  arranging  it  in  quadratic  form.  The  first  term 
in  Eq.  197  is  real,  and  represents  the  normal  mode  of  the  system  in  the  absence  of  driving 
forces.  The  quantities  under  the  square  root  sign  determine  the  stability  of  the  system.  The 
first  term  under  the  square  root  sign  represents  the  difference  in  the  centrifugal  force  term 
across  the  interface  and  surface  tension  between  the  liquid  metal  and  cover  gas.  Both  terms 
are  stabilizing  for  the  case  of  a  liquid  metal  supported  by  a  less  dense  cover  gas.  The 
second  term  under  the  square  root  is  always  destabilizing  and  is  the  source  of  the  Kelvin- 
Helmholtz  instability  and  the  discontinuity  in  primary  flow  across  the  interface. 

In  Eq.  197,  if  the  quantity  Q  under  the  square  root  [Q]'A  is  greater  than  0,  Q  >  0, 
the  angular  frequency  a>  is  real  and  neutral  stability  is  maintained  in  the  system.  However, 
if  Q  <  0,  then  c o  =  a  ±  ib  and  the  fluid  flow  system  becomes  unstable.  In  general, 
because  of  the  dependence  of  quantities  on  the  wave  propagation  constant  &  in  Eq.  197, 
instabilities  occur  in  the  region  of  the  k  spectrum  when  the  sum  of  the  quantities  under  the 
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square  root  sign  is  negative,  i.e.,  Q  <  0.  Thus,  the  stability  threshold  is  obtained  from 


( 


QSl  QgVg  +  k2T\  >  {Vg  _  Vl)2  t 


R  R 

which  reduces  to 

<Yg  ~  VL) 2  <  — 
Qc 


(Gg+Gz.) 
kT 


V  M  Ql  ) 


(198) 


(199) 


when 


Gg^C  : 


is  neglected  with  respect  to 


QiYl 


and  set  qg  +  ql~  ql. 


The  wave  number  for  the  first  unstable  mode  for  the  system  k  is  derived  by  finding 
the  minimum  value  of  the  right  hand  side  of  Eq.  199: 


For  this  k,  the  stability  threshold  equation  (Eq.  199)  is  expressed  as 

(VG  ~  VL)2  <  —  •  (200) 

Go*  * 

If  the  azimuthal  velocity  profiles  similar  to  those  in  Fig.  6  are  assumed  for  the  primary 
flow,  the  following  approximations  can  be  made  where  Vow  is  the  rotor  velocity  at  the 
interface: 


(VG  -  VL)2  =  a2V0\  , 

(201) 

and 

vL  =  - 

(202) 

where  a2  =  (1  . 

i  *  AXIAL  DISTANCE 


Fig.  6.  Assumed  typical  azimuthal  flow  profile  Va(z)  at  interface  (r  =  R)  based  on  similarity  solution 
in  base  flow  and  conservation  of  angular  momentum  at  the  interface.1 
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Using  these  approximations  in  Eq.  200,  the  critical  speed  for  ejection  is  obtained  as 


^CRITICAL  ~  (2/?>CRJTICAL 


2  (QLT) 


'h 


and 


‘CRITICAL 


2(QLT) 


Qg^ 


3/2 


(203) 


(204) 


(/critical  is  the  rotor  speed  at  the  interface  in  meters  per  second,  and  ^critical  is  the 
critical  angular  velocity  at  which  hydrodynamic  instability  takes  place.  These  results  agree 
exactly  with  those  obtained  by  Woo  et  a/.,1  by  a  somewhat  different  derivation. 
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